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Abstract 

Decay constants of P-wave mesons are computed in the framework of instantaneous Bethe- 
Salpeter method (Salpeter method). By analyzing the parity and possible charge conjugation 
parity, we give the relativistic configurations of wave functions with definite parity and possible 
charge conjugation parity. With these wave functions as input, the full Salpeter equations for 
different P-wave states are solved, and the mass spectra as well as the numerical values of wave 
functions are obtained. Finally we compute the leptonic decay constants of heavy-heavy and 
heavy-light ^Pq, ^Pi and ^Pi states. 
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I. INTRODUCTION 



The quantities of meson decay constants are important, since they play important roles 



in many aspects, and the studies of them have become hot topics in recent years 



aaay 
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12l |. But most of these investigations are focused on the estimating decay 



edge of decay constants for P-wave 
|l3|,Q,[l5|,Q,H Il8l | . We present a careful 



constants for S-wave mesons, and we lack the know^ 

study of the decay constants for heavy ^Pq, ^Pi and ^Pi states including the relativistic 
corrections. 

In previous letters 19| , the decay constants of heavy- heavy and heavy-light pseudoscalar 
(^^o) and vector {^Si) mesons are calculated in the framework of relativistic instantaneous 
Bethe-Salpeter method [20] (also called Salpeter method [21]), good agreement of our pre- 
dictions with recent lattice, QCD sum rule, other relativistic model calculations as well as 
available experimental data is found. 

In this letter, we extend our previous analysis to include P-wave mesons, present the 
calculations of decay constants for heavy P-wave states in the framework of full Salpeter 
equation which is a relativistic method. Based on the S — L coupling scheme, we analyze the 
parity and possible charge conjugation of ^Pq, ^Pi and ^Pi bound states, and give general 
formula for the wave functions which are in relativistic form with definite parity and charge 
conjugation symmetry (0"*'"'', 1^^, 1^* for equal mass ^Pq, ^Pi and ^Pi bound states, and O'^, 
l'^, 1+ for un-equal mass bound states). Then with these wave functions as input, we solve 
the full Salpeter equations, obtain the mass spectra and numerical values of wave functions 
for different P-wave states. Finally, we compute the leptonic decay constants for heavy ^Pq, 
^Pi and ^Pi states. 

This letter is organized as following, in section II, we introduce the relativistic Bethe- 
Salpeter equation and Salpeter equation. In section III, we give the formula of relativistic 
wave functions and decay constants for P-wave states. We solve the full Salpeter equations, 
obtain the mass spectra and wave functions of P-wave mesons. Finally, we use these rela- 
tivistic wave functions to calculate the decay constants of heavy P-wave mesons and show 
the numerical results and discussion in section IV. 
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II. INSTANTANEOUS BETHE-SALPETER EQUATION 



In this section, we briefly review the instantaneous Bethe-Sal 



Deter equation and introduce 
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our notations, interested reader can find the details in Ref. 
TheBethe-Salpete.(BS)e,ua«on..ead.Q: 

(^1 - mi)x(g)(j^2 + m2)=ij J^^y^P, q)x{k) , (1) 

where xi^l) is the BS wave function, V{P, k, q) is the interaction kernel between the quark 
and antiquark, and pi,p2 are the momenta of the quark 1 and anti-quark 2. The total 
momentum P and the relative momentum q are defined as: 

Pi = aiF + g, «! — 



P2 = - q, a2 



nil + ^2 
m2 



nil + "^2 

We divide the relative momentum q into two parts, q\\ and q±, 

q' = qi; + ql, 

gf ^ (P ■ q/M')P^ , ^ gA< - gf . 

Correspondingly, we have two Lorentz invariant variables: 



(ip = ^ ^ ^T = \Jql-q^ = \f-qi ■ 

When P= 0, they turn to the usual component go and |g|, respectively. 

After instantaneous approach, the kernel V{P^ k, q) takes the simple form: 

V{PXq)^V{k^,q^). 

Let us introduce the notations (Pp{qj_) and ri{qj_) for three dimensional wave function as 
follows: 

dqp^ 



Viql)^ J ^,V{k^,qi.)¥>p{kl) . (2) 



(2vr) 

Then the BS equation can be rewritten as 



xiq\\,q±) = Siipi)7]iq±)S2iP2) ■ (3) 
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The propagators of the two constituents can be decomposed as: 



Si{Pi) 



J{i)qp + aiM - Ui + ie J{i)qp + a-M + Ui - ie ' 



with 



2u. 



L 

M 



uJi ± J{i){mi + 



(4) 



(5) 



where i = 1,2 for quark and anti-quark, respectively, and J{i) = (—1)*^^. Here A.fp{q±) 
satisfy the relations: 

K%{q±) + Arp(g±) = ^ , A±(g^)^Aj(g^) = K%{q^) , Aj(g^)^A;(g^) = . (6) 



Introducing the notations as: 



(7) 



and we have 

With contour integration over qp on both sides of Eq.(l3]), we obtain: 

. A^p(gx)r/p(gx)A^p(gx) Arp(g±)r7p(g±)A2p(g±) 



and the full Salpeter equation: 



(M — a;i — (^2) 



(M + + UJ2) 



{M~uJi- uj2)^++{q^) = A+ (gx)r?p(g±)A+ (gx) , 
{M + UJI+ uj2)(fp~iq±) = -Aipiq±)Vp{q±)^2pi(l±) 

¥^p~{q±) = fp^{g±) = . 

The normalization condition for BS wave function is: 



27r2 



2Pn . 



(8) 



(9) 



In our model, Cornell potential, a linear scalar interaction plus a vector interaction is 
chosen as the instantaneous interaction kernel V 221 ]: 

n^) = K(g1+7o®7°K(g1, 

1 2 asiq) 



Vsiq) = -i- + Vo)6'iq) + -^ . , 

where the coupling constant as{q) is running: 

127r 1 



37r2 (g^ + «2) 



(10) 



(^s{q) 



27 log(a + -^ 





^QCD' 



and the constants A, a, a, Vq and Aqcd are the parameters that characterize the potential. 



III. RELATIVISTIC WAVE FUNCTIONS AND DECAY CONSTANTS 



In this section, by analyzing the parity and possible charge conjugation parity of corre- 
sponding bound state, we give a formula for the wave function that is in a relativistic form 
with definite parity and possible charge conjugation parity symmetry. 

A. Wave function for ^Pq state 

The general form for the relativistic Salpeter wave function of ^Pq state, which = O"*" 
(or J-^*^ = 0''''*' for equal mass system), can be written as: 



^oAq±) = fi{q±)^+ /2(?±)^ + + h{q±)r- (n) 

The equations 

cp+-{q^) = ^-+{q^) = (12) 
give the constraints on the components of the wave function, so we have the relations 

j^^^^-^ ^ /i(g±)gi(mi + m2) ^ /2(9±)gi(^i -UJ2) 

M{ujiuj2 + mim2 + gi) ' M{miLJ2 + m2UJi) 

Then there are only two independent wave functions fi{q±) and f2{(l±) been left, from Eq. 
(IH]), we obtain two coupled integral equations, by solving them we obtain the numerical 
results of mass spectra and wave functions, interesting reader can find the details of this 
method in Ref. [2^ or Ref. 23]. 



In our calculation, we choose the center-of-mass system of the corresponding state, so 
q\\ and q± turn to the usual components {qo,0) and (0,g), cui = (m^ + cfY^'^ and UJ2 = 
{m\ + g^)^/^. The normalization condition for the ^Pq wave function is: 

^ dg^^l6A^W^2M. (13) 
B. Decay constant of ^Pq state 

The decay constant Fsp^, of scalar ^Pq meson is defined as 

(0|gi7,.(l - 75)g2|'Po) = Fsp,P„ (14) 
which can be written in the language of the Salpeter wave functions as: 



5 



(0|g-i7M(l - 75)?2|'Po) ^^cjTr [(/^o+(?)7m(1 - Ts)] 



(27r) 
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Therefore, we have 



^ ^ ^ M a;iC(;2 + m\m2 — (p M{miUJ2 + ?«2'^i) / ^'^ 

= 4^/iv/^(' (15) 

^V^'^y (27r)3 1, M(miu;2 + m2u;i)j ^ ' 



4v^ /■ /2g"(a;2-c^i) .^g. 
7 f27r)3 fmia;2 + moo^i) 



° M J (27r)3 (miu;2 + m2u;i) 



C. Wave function for ^Pi state 

The general form for the Salpeter wave function of ^Pi state, which — 1+ (or J^'-^ — 
1++ for equal mass system), can be written as: 

vMqi.) = ts^^^pP^qle" [hM^ + 12^1^ + h4i_l^ + lUe^""' qLpPalsl,/ m] /M\ (17) 
The equations 

vt^{(lA.) = V-^^iqi.) = (18) 

give the constraints on the components of the wave function 

r, . fi{q±)M{miU2 - m2UJi) . f2{q±)M{-uiUJ2 + mim2 + qj_) 

J3[q±) = 27 ] \ ' /4l?-Lj = 27 ] \ 

qi{uji + u;2) ql{mi + m2) 

The normalization condition for the ^Pi wave function is: 

dq 32/i/2CJicj2(c^ic^2 - "^i"^2 + g^) ^ 2_/vf (19) 



(27r)3 3(mi + m2)(co'i + Ci;2) 

D. Decay constant of state 

The decay constant Fsp^ is defined as 

(01^17^(1 -75)g2|'Pi,e) = Fap^MeJ, (20) 
and can be formulated using the Salpeter wave function as: 

(01^17^(1 -75)g2|'Pi, 6) = ^|Ncj rr[(^i+(a7^(l-75)]^, (21) 

then we have 



(27r)3' 

^ 8A/iVe /■ f2{uJlUJ2 - 17111712 + g^) . s 

~ 3M J (27r)3 (mi + mo) ' ^ ^ 
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E. Wave function for ^Pi state 



The general form for the Salpeter wave function of ^Pi state, which = l"*" (or J^^ 
1+^ for equal mass system), can be written as: 



<^i+(g±) = g± ■ el 



p J 

fl{q±) + + ^3(9±)^ + /4(9±y^2 



75- (23) 



The equations 

¥^^-(g±) = ¥^r+^(g±) = (24) 

give the constraints on the components of the wave function 

rr ^ /i (g± )M(mi - ma) , f2{q±)M{uJi + U2) 

J3[q±> = --, [ 2T' /4WJ - 



{uJiU2 + mim2 - q\) ' (micua + macui 

The normalization condition for the ^Pi wave function is: 

dq 16fif2UJiUJ2^ 



/ / = 2M. (25) 

(27r)-^ 3{miUJ2 + m2UJi) 

F. Decay constant of ^Pi state 

The decay constant Fip^ is defined as 

(0|gi7M(l-75)g2rPi,e) = Pp.Me^^, (26) 
which can be formulated as: 



(0|gi7M(l-75)g2rPi,e) = / Tr (5)7^(1 - 75)] (27) 
finally, we obtain 



^ iy/K f dq h{mi - m2)q^ , . 

3M J (27r)3 + mims + g^) ■ ^ ^ 

IV. NUMERICAL RESULTS AND DISCUSSION 

In our method, there are some input parameters appearing in the potential, we need to fix 
them when solving the full Salpeter equations. Usually, we fixed the parameters by fitting 



the experimental mass spectra for mesons, but for P-wave states, we 
so we adopt almost the same parameters as in the 0~ states Ref. 



ack experimental data. 



22], and only vary the 



parameter Vq by fitting the ground P-wave cc states, Xco, Xci and he- In previous letter [l9|, 
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we found if we choose same parameters set, the mass predictions of our model can not agree 
very well with experimental data for pseudoscalar and vector mesons, we find the same thing 
happens to the different P-wave states, so we vary the only possible different parameter Vq 
to fit the data. For ^Pq states, we choose Vq = —0.566 GeV, for ^Pi states, Vq = —0.452 
GeV, and for ^Pi states, Vq = —0.437 GeV. The values of other parameters are same as in 



the states Ref. 
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a = e = 2.7183, a = 0.06GeV, A = 0.20GeV^ Aqcd = 0.26GeV and 
rub = 5.224GeV,mc = 1.7553GeV,m, = 0.487GeV, = 0.311GeV,m„ = 0.305GeV. (29) 



TABLE I: Mass spectra in unit of MeV for cc and bb P wave states. 





Ex(cc) 


cc 


Ex{bb) 


bb 


1 '^'o 


3415.2 


3415.9 


9859.9 


9860.1 


2 3po 




3831.1 


10232.1 


10223.9 


3 3Po 




4132.4 




10497.0 


4 3po 




4369.4 




10719.2 


1 ^Pi 


3510.6 


3510.9 


9892.7 


9892.1 


2 




3923.1 


10255.2 


10255.0 


3 




4222.0 




10527.4 


4 




4456.9 




10750.0 


I 'Pi 


3524.4 


3524.4 




9900.4 


2 




3935.8 




10262.6 


3 




4234.2 




10534.6 


4 




4468.7 




10757.1 



We show our theoretical predictions of mass spectra for cc states up to the 4P states as 
well as the experimental data in Table I. One can see that, our predictions for mass splitting 
are, 2P - IP ~ 410 MeV, 3P - 2P ~ 300 MeV, 4P - 3P ^ 235 MeV; the predicted mass 
for state 2^Pi is 3923 MeV, this is a little smaller but consist with the traditional prediction 
of potential model , which is about 50 MeV higher than the one of new state X(3872). We 
show the predicted mass spectra for other states In Table II, the interesting quantity is also 
the mass splitting between the first radial excited state and ground state, for all the ^Pq, 
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and ^Pi states, 2P - IP ~ 330 MeV for ub, 2P - IP ~ 345 MeV for uc, 2P - IP ^ 362 
MeV for sb and 2P - IP ~ 381 MeV for sc. 



TABLE II: Mass spectra in unit of MeV for heavy P wave states. 





1 3 D 

t Pq 


O 3 ZD 

2 Pq 


1 3 D 




T 1 P 
1 Pi 


2 Pi 


cb 


6728.7 


7127.8 


6829.5 


7225.3 


6845.1 


7239.6 


sb 


5767.2 


6130.3 


5830.9 


6192.2 


5836.4 


6197.3 


db 


5667.9 


5998.9 


5711.7 


6042.8 


5709.2 


6041.4 


ub 


5664.8 


5994.1 


5707.6 


6037.2 


5704.7 


6035.5 


sc. 


2386.5 


2767.4 


2447.8 


2827.3 


2449.8 


2830.4 


dc 


2273.1 


2619.2 


2314.1 


2661.3 


2307.6 


2657.7 


uc 


2269.3 


2613.7 


2309.5 


2655.0 


2302.5 


2651.1 



We also calculate the mass spectra for P-wave bb system, as argued in Ref. [2J], there 
are double heavy b quarks, and the flavor Nf = 4, so we have to choose a new set of 
parameters as well as smaller value of coupling constant. We change the previous sea 



parameters to Aqcd = 0.20 GeV, nib = 5.13 GeV which have been adopted in Ref. |24| . 
choose Vo = -0.553 GeV for ^Pq states, Vq = -0.521 GeV for ^Pi states, Vo = -0.514 
GeV for ^Pi states, and other parameters are not changed. With this set of parameters, 
the coupling constant at the scale of bottom quark mass is as{mb) = 0.23. The numerical 
results and experimental data of mass spectra for bb system are also shown in Table I. One 
can see that our predictions, 2P — IP ^ 363 MeV, can fit the experimental data very well, 
and our mass splitting prediction, 3P - 2P ~ 273 MeV, 4P - 3P ~ 223 MeV. 

TABLE III: Decay constants in unit of MeV for P wave cc and bb states. 





n 


n 


1 


2 ^Pi 


3 


4 ^Pi 


cc 








206 


-207 


199 


-189 


bb 








129 


-131 


126 


-121 



Becides the mass spectra, we also obtained the relativistic wave functions for heavy 
mesons when solving the full Salpeter equation. With these wave functions, we calculated 
the decay constants for heavy-heavy and heavy-light P wave mesons. In Table III, we 
show our estimates of decay constants for cc and bb ^Pi-wave systems up to third radial 
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TABLE IV: Decay constants in unit of MeV for heavy P wave states. 



1 2 1 ^p^ 2 ^Pi 1 2 

d) 88 -85 160 -165 50 -49 

sb 140 -130 157 -156 76 -71 

d5 145 -129 150 -144 76 -70 

ub 145 -128 150 -143 76 -70 

sc 112 -91 219 -204 62 -50 

dc 132 -102 212 -190 72 -56 

uc 133 -102 211 -189 72 -56 



excited states, for ^Pq and ^Pi equal-mass states, the decay constants vanish. The relative 
sign between the Fip and F2P are minus, this come from the behavior of wave function, 
since there is a node in the 2P radial wave function, the minus sign means the dominant 
contribution come from the part after the node. Our predictions show that, the decrease of 
the numerical value of decay constant for higher excited state is not evident comparing with 
the lower excited state, for example, F^p — Fip = 7 MeV for ^Pi cc system, and F^p—Fip = 3 
MeV for ^Pi bb system. 



TABLE V: Decay constants in unit of MeV for heavy P wave states in different models. 







1 






1 






1 






ours 


[13] 


[14] 


ours 


[13] 


[14] 


ours 


[13] 


fl4] 


sb 


140 




146 


157 




181 


76 




84 


ub 


145 


112 


162 


150 


123 


187 


76 


68 


93 


sc 


112 


71 


110 


219 


121 


240 


62 


38 


63 


uc 


133 


86 


139 


211 


127 


249 


72 


45 


82 



In Table IV, we show our estimates of decay constants for unequal-mass P-wave ground 
and first radial excited states. It is observed that the decay constants of ^Pi states are much 
larger than those of the corresponding ^Pi states, while the corresponding values for ^Pq is 
between them. 

For comparison, we show our predictions for decay constants and other theoretical pre- 
dictions jisl . [1^ in Table V. We have changed results in Ref. [l^ from the j — j coupling 
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scheme to S — L coupling scheme by using the following equations |13l. |25| 



= ipf/^) - ^ \pn i^Pi) = ^ lA^/^) + ^jl \pn (30) 

Rough agreement can be found between the values of decay constants estimated by different 
methods, this means we need more effort for the knowledge of P-wave decay constants. 

In conclusion, we estimated the decay constants for heavy P-wave ^Pq, ^Pi and ^Pi 
mesons in the framework of the relativistic Bether-Salpeter method. 

This work was supported in part by the National Natural Science Foundation of China 
(NSFC) under Grant No. 10675038. 
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